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We have measured noise in thin-film superconducting coplanar waveguide resonators. This noise 
appears entirely as phase noise, equivalent to a jitter of the resonance frequency. In contrast, 
amplitude fluctuations are not observed at the sensitivity of our measurement. The ratio between 
the noise power in the phase and amplitude directions is large, in excess of 30 dB. These results 
have important implications for resonant readouts of various devices such as detectors, amplifiers, 
and qubits. We suggest that the phase noise is due to two-level systems in dielectric materials. 
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Thin-film superconducting microwave resonators are of 
interest for a number of applications, including the multi- 
plexed readout of single electron transistors (SET)|l|, mi- 
crowave kinetic inductance detectors (MKID)^,^ nor- 
mal metal-insulator-superconductor (NIS) tunnel junc- 
tion detectors superconducting quantum interference 
devices (SQUID) [1,11], and qubits 0, The device to 
be measured presents a variable dissipative or reactive 
load to the resonator, influencing the resonator qual- 
ity factor Qr or frequency fr, respectively. Changes to 
both Qj. and fr may be determined simultaneously by 
sensing the amplitude and phase of a microwave probe 
signal[2]. While several early demonstrations used hand- 
assembled lumped-element circuits P, |3, Q, frequency- 
domain multiplexing of large arrays generally will re- 
quire compact microlithographed high-Q^ resonators J|. 
Such resonators are also needed for strong coupling to 
charge qubits [?!]. Noise in microlithographed resonators 
has been observed [2, ^ and can be a limiting factor for 
device performance, but is not well understood. In this 
letter, we report measurements of resonator noise, show 
how the noise spectra separate into amplitude and phase 
components, and discuss the physical origin of the noise. 

We studied quarter-wavelength coplanar waveguide 
(CPW) resonators^ (Fig. la) with center strip widths 
w of 0.6 to 6 /im and gaps g between the center strip 
and ground planes of 0.4 to 4 /xm, and with impedances 
Zo « 50 ri. Resonator lengths of 3 to 7 mm produce res- 
onance frequencies fr between 4 and 10 GHz. Frequency 
multiplexed arrays of up to 100 resonators are coupled 
to a single CPW feedline. The CPW circuits are pat- 
terned from a film of either Al (T^ = 1.2 K) or Nb (T^ 
= 9.2 K) on a crystalline substrate, either sapphire. Si 
or Ge. The surfaces of the semiconductor substrates are 
not intentionally oxidized, although a native oxide due 
to air exposure is expected to be present. 

A microwave synthesizer at frequency / is used to ex- 
cite a resonator. The transmitted signal is amplified with 
a cryogenic high electron mobility transistor (HEMT) 
amplifier and is compared to the original signal using 
an IQ mixer, whose output voltages / and Q are propor- 
tional to the in-phase and quadrature amplitudes of the 
transmitted signaljl, 3 (see Fig. 2 inset). As / IS var- 
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FIG. 1: (Color online) (a) Schematic illustration (not to scale) 
of the resonator and feedline geometry; see Day et aZ. 0] for the 
equivalent circuit. Black represents the superconducting film; 
white represents bare substrate. The coupler relies on the mu- 
tual capacitance between the center strips of the feedline and 
resonator CPW lines and may be considered to be a lumped 
element C = 1/ (/r^o VSt^Qc) where Qc is the coupling- 
limited quality factor, (b) Resonance circle of a 200 nm Nb on 
Si resonator at 120 mK (solid line), quasiparticle trajectory 
calculated from Mattis-Bardeen theory [13] (dashed line). For 
this figure, the readout point ^ = ]/, Q] is located at the res- 
onance frequency fr. (c) Noise ellipse (magnified by a factor 
of 30). Other parameters are /^=4.35 GHz, Qr = 3.5 x 10^ 
(coupling limited), it)=5 ^m, 17=1 /^m, readout power P,, ~ 
-84 dBm and internal «-30 dBm. 



ied, the output ^ = ]/, QY" (the superscript T represents 
the transpose) traces out a resonance circle (Fig. 1(b)). 
With / fixed, ^ is seen to fluctuate about its mean, and 
the fluctuations 5S,{t) — [5I{t),6Q{t)]'^ are digitized for 
noise analysis, typically over a 10 s interval using a sam- 
ple rate of 250 kHz. 

The fluctuations 5S,{t) are observed to be primarily in 
the direction tangent to the resonance circle, while the 
fluctuations in the orthogonal direction are small. These 
two directions correspond to fluctuations in the phase 
and amplitude of the resonator's electric field E, respec- 
tively. This observation can be quantified by studying 
the spectral-domain noise covariance matrix S(y), de- 
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FIG. 2: Noise spectra in the phase {Saa{v), sohd hne) and 
amphtude (Sbbii^), dashed Une) directions, and the rotation 
angle {(piv), dotted hne). The noise data are from the same 
Nb/Si resonator under the same condition as in Fig. 1. The 
inset shows the diagram of the homodyne readout system. 



fined by 

(1) 

where S^i^i') is the Fourier transform of the time-domain 
data, the dagger represents the Hcrmitian conjugate, 
Sii{v) and Sqq{v) are the auto-power spectra, and 
^igii^) is the cross-power spectrum. The matrix Siiy) 
is Hermitian and may be diagonahzed using a unitary 
transformation; however, we find that the imaginary part 
of SiQ is neghgible and that an ordinary rotation apphed 
to the real part Re 5(1/) gives almost identical results. 
The eigenvectors and eigenvalues are calculated at every 
frequency, i^: 
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where 0{i') = [vaiv), Vb{v)\ is an orthogonal rotation ma- 
trix. We use Saa(y) and Va(v) to denote the larger eigen- 
value and its eigenvector. 

A typical pair of spectra Saa{v) and Shb{v) are shown 
in Fig. 2, along with the rotation angle (t){v), defined as 
the angle between Va{v) and the / axis. Three remarkable 
features are found for all noise data. First, ^(i/) is inde- 
pendent of V within the resonator bandwidth (the r.m.s. 
scatter is < 0.4° per 10 Hz frequency bin from 1 Hz 
to 5 kHz), which means that only two special directions, 
Va and Ub, diagonalize S{iy). Eqn. ^ shows that Saa{i^) 
and Sbb{t^) are the noise spectra projected into these two 
constant directions. Second, Va is always tangent to the 
IQ resonance circle while Vb is always normal to the cir- 
cle, even when / is detuned from fr- Third, Saai^^) is 
well above Sbbii^) (see Fig. 2). The character of the noise 
can be clearly visualized by plotting a noise ellipse, de- 
fined by = 1, where C = Q ReS(y)dv is the 
covariance matrix for 51 and 5Q filtered for the corre- 
sponding bandpass (we use vi = 1 Hz and 1^2 = 1 kHz). 
The major axis of the noise ellipse is always in the phase 



direction, and the ratio of the two axes is always very 
large (Fig. 1(c)). 

Fig. 2 also shows that the amplitude noise spectrum is 
flat except for a, l/u knee at low frequency contributed 
by the electronics. The amplitude noise is independent 
of whether the synthesizer is tuned on or off the reso- 
nance, and is consistent with the noise temperature of 
the HEMT amplifier. The phase noise spectrum Q has 
a.l/v slope below 10 Hz, typically a v^^l"^ slope above 
10 Hz, and a roll-off at the resonator bandwidth fr/2Qr. 
The phase noise is well above the HEMT noise, usually 
by two or three orders of magnitude (in rad^/Hz) at low 
frequencies. It is well in excess of the synthesizerphase 
noise contribution or the readout system noise [l,^. 

Quasiparticle fluctuations in the superconductor can 
be securely ruled out as the source of the measured noise 
by considering the direction in the IQ plane that would 
correspond to a change in quasiparticle density Suqp. 
Both the real and inductive parts of the complex con- 
ductivity a respond linearly to Suqp, Sa — 5ai — iSa2, 
resulting in a trajectory that is always at a non-zero an- 
,gle ip = tan~^(5cri/(5(T2) to the resonance circle, as indi- 
cated by the dashed lines in Fig. 1(b) and (c). Mattis- 
BardeenH^ calculations yield ip > 7° for Nb below 1 K, 
so quasiparticle fluctuations are strongly excluded since 
-0 >> (70. Furthermore, "0 is measured experimentally 
by examining the response to X-ray, optical/UV, or sub- 
millimeter photons, and is typically tp ss 15°[T]|. 

The phase noise depends on the microwave power in- 
side the resonator (Pint), the materials used for the res- 
onator, and the operating temperature. The power de- 
pendence for various material combinations is shown in 
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Fig. 3; all follow the scahng Saa{v) oc Pj„j . For compar- 
ison, amplifier phase noise is a multiplicative effect that 
would give a constant noise level independent of Pinti 
while the amplifier noise temperature is an additive ef- 
fect that would produce a 1/Pint dependence. Sapphire 
substrates generally give lower phase noise than Si or Ge. 
However the Nb/Si device showed low noise comparable 
with Al/sapphire, suggesting that the etching or inter- 
face chemistry, which is different for Nb and Al, may 
play a role. Two Al/Si resonators with very different Al 
thicknesses and kinetic inductance fractions [l2] fall onto 
the dashed equal-noise scaling line, strongly suggesting 
that the superconductor is not responsible for the phase 
noisef^. Furthermore, the noise of a Nb/Si resonator de- 
creased by a factor of 10 when warmed from 0.2 K to 1 K, 
more strong evidence against superconductor noise since 
Nb has Tc — 9.2 K and its properties change very little for 
T << Tc- More detail on the temperature dependence 
will be published separately. 

The evidence leads us to suggest that the noise is 
caused by fluctuating two-level systems (TLS) in the 
dielectric materials — either the bulk substrate or its 
exposed surface, the interface layers between the metal 
films and the substrate, or any oxide layers on the metal 
surfaces. Models assuming a collection of TLS with a 
wide range of excitation energies E and relaxation rates 
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FIG. 3: Power and material dependence of the phase noise 
a.t V = 1 kHz. To compare resonators with different fr and 
Qr, phase noise is converted to fractional frequency noise, 
calculated by Ssj\,{v) / = Saa{y) /'^Ol- AH the resonators 
have w=3 /xm, g=2 fj,m and are measured around 120 mK. 



have long been used to explain the low t emp erature phys- 
ical properties of non-crystalline solids flSl ll^ ITsI . TLS 



are also found in crystalline materials j id 17 1 but at 
lower densities. Fluctuations due to TLS are of inter- 
est theoretically [H, [l^ and have been observed in di- 
electrics, either as telegraph or l/v noise, using tunnel 
junctions [13], sin gle electron transistors 21 1 , and atomic 
force microscopes [22|. Other examples include telegraph 
noise in the resistance of metallic nanoconstrictions |23j 
and qubit dephasing effects Q. 

Recent experiments 0, have illustrated the reactive 
loading effect that occurs when a single TLS (a qubit) 
is coupled to a microwave resonator. For weak coupling 
g or large detuning A/ = \E/h — fr\, g « A/, this 
reactive loading causes the resonator frequency to shift 
by ±g^/A/ depending on the quantum state of the TLS. 
Thermal fluctuations - absorption or emission of thermal 
phonons by a collection of TLS - could therefore cause 
phase noise, equivalent to a fluctuating dielectric con- 
stant 5e{f, t) . For this model one expects the noise to van- 
ish exponentially for T << hfr/2k as the TLS settle into 
their ground states, while the increase in TLS-phonon 
transition rates might explain the observed noise decrease 
at high temperatures. However, at present we do not 
have a model which quantitatively explains the data. 
Alternatively, TLS frequency fluctuations SE{t)/h pro- 
duced by TLS-TLS interactions [isl [25j and observed in 
single molecule optical fluorescence experiments [1^ [23| 
should give phase noise with a power-law rather than 
exponential decrease at low temperatures. 
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The Pint noise scaling is indicative of TLS satura- 
tion; otherwise, the phase noise would be independent 
of Pint as expected for dielectric constant fluctuations. 
Similarly, power-independent thermal fluctuations of the 
TLS dielectric polarization SP [19| are ruled out since 



ri^Ti is independent of Aq. 



this would give additive noise with equal amplitude and 
phase components and scaling as 1/Pint- TLS satura- 
tion effects are well known [sT [l5l|: indeed, fr and 
show anomalous temperature and power dependence [28| 
for T << Tc- Saturation of a single TLS depends on 
its frequency detuning A/, its position r in the spa- 
tially varying resonator field E(f), and the orientation 
and strength of its dipole moment d. The growth of the 
saturation zone in this parameter space with increasing 
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Pint may explain the observed Pj„t noise scaling. 

The microwave fields in our resonators are sufficiently 
strong to cause TLS saturation. In the Bloch equa- 
tion framework, saturation for zero detuning occurs when 
^^TiT2 > 1, where = d ■ E/h is the Rabi frequency 
and Tj~i, T^^ are the usual energy relaxation and de- 
phasingrates[29| . The distribution of T^^ is extremely 
broad [15[, extending above and below the range of noise 
frequencies we observe, because T-^^ oc Aq/E^ is con- 
trolled by the tunnel coupling Aq and is therefore expo- 
nentially sensitive to the potential barrier (Tsl [29j . For 
silica, Tf 1 « 1 MHz(£;/4GHz)3 {Aq/ E)"^ coth{E/2kT). 

The dipole moment d is proportional to Ao/-E[iB|, so 

Meanwhile, scales as 
and has a relatively narrow distribu- 
tion; for silica, T^' « lOMHz (T/200mK)i-5[305. Satu- 
ration is therefore mostly independent of Ag, and occurs 
when E ■ d/\d\ exceeds a critical value Edfr, T). Indeed, 
the shape of the observed phase noise spectrum does not 
change significantly with Pint until the superconductors 
become nonlinear^. For silica, i?c(4 GHz, 200 mK) « 
31Vm~i and corresponds to Pint ~ — 75 dBm which is 
significantly below the power levels we use (see Fig. 3). 
We also find Ec{7.2 GHz, 25 mK) « lOVm-i, which is 
appropriate for the experimental conditions of Martinis 
et aL@| and is similar to their observed onset of TLS 
saturation. 

In summary, we find that lithographed resonators dis- 
play v~^l'^ noise directed purely in the phase direction 
which varies with readout power, temperature, and the 
substrate material. These results are important for the 
device optimization — devices relying on resistive load- 
ing should be able to avoid this noise source. For MKID 
detectors, the sin'0 amplitude component of the signal 
is already available when using an IQ readout, and an 
optimally weighted phase/amplitude measurement can 
substantially improve the sensitivity at low frequencies 
where the phase noise is larger. 
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